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SUMMARY 


Equilibrium convective heat transfer in several real gases was investi- 
gated. The gases considered were air, nitrogen, hydrogen, carbon dioxide, and 
argon. Solutions to the similar form of the boundary -layer equations were 
obtained for flight velocities to 30,000 ft/sec for a range of parameters 
sufficient to define the effects of pressure level, pressure gradient, 
boundary -layer -edge velocity, and wall temperature. Results are presented 
for stagnation -point heating and for the heating -rate distribution. 

For the range of parameters investigated the wall heat transfer depended 
on the transport properties near the wall and precise evaluation of properties 
in the high-energy portions of the boundary layer was not needed. A corre- 
lation of the solutions to the boundary-layer equations was obtained which 
depended only on the low temperature properties of the gases. This result can 
be used to evaluate the heat transfer in gases other than those considered. 



The largest stagnation -point heat transfer at a constant flight velocity 
was obtained for argon followed successively by carbon dioxide, air, nitrogen, 
and hydrogen. The blunt-body heating-rate distribution was found to depend 
mainly on the inviscid flow field. 

For each gas, correlation equations of boundary-layer thermodynamic and 
transport properties as a function of enthalpy are given for a wide range of 
pressures to a maximum enthalpy of 18,000 Btu/lb. 

INTRODUCTION 



Convective heat transfer to bodies entering the earth f s atmosphere has 
been studied extensively. No attempt will be made to reference all these 
investigations, but representative examples are references 1 to 3. Only a few 
investigators have examined the problem associated with flight into planetary 
atmospheres having constituent gases differing from air and these deal 
directly with stagnation -point heat transfer. (See, e.g., refs. 4 and 5*) 

Prior to these latter investigations, it was common to correlate heat -transfer 
results in terms of the transport and thermodynamic properties at the boundary 
layer edge (i.e., see ref. 2 ), but as indicated in reference 5; a correlation 
for various dissociating gases was obtained in terms of molecular weight of 
the cold mixture. A correlation of this nature is much more convenient since 
estimates of heating in any number of gases seemed to be feasible without 
requiring the burdensome task of evaluating the transport properties at the 
boundary-layer edge. In light of these two apparently different modes of cor- 
relation, it seemed appropriate to study the relation of the transport proper- 
ties to the heat transfer. 


In a cursory look at this problem, pertinent thermodynamic and transport 
properties of several gases were compared because it is through these proper- 
ties that differences in heat transfer would be expected to appear. This com- 
parison illustrated some of the property differences between gases and these 
differences were investigated further to assess their effects on the convec- 
tive heat transfer. 

It is the purpose of this report to present the results of this investi- 
gation for several real gases, including air, and to point out the significant 
differences and similarities between the results for the various gases. The 
results were obtained by solving the boundary -layer equations for the gases, 
air, nitrogen, hydrogen, carbon dioxide, and argon, subject to the assumptions 
of local similarity and thermodynamic equilibrium. Solutions are presented 
for flight velocities to 30,000 ft/sec and for a range of parameters suffi- 
cient to define the influence of pressure, pressure gradient, wall temperature, 
and velocity at the edge of the boundary layer. The results are presented for 
stagnation -point heating and heating -rate distribution and are correlated in 
terms of the gas properties at low temperatures. 


SYMBOLS 


C P 


c 


Pi 


Ci 


D 


ij 


f 

g 

h 

H 

H s 

k 

k f 


n 

n 


total specific heat 
specific heat of species i 
mass fraction of species i 
multicomponent diffusion coefficient 


'L 


u , 


normalized total enthalpy, 
static enthalpy 


JL 

He 


total enthalpy 



stagnation enthalpy 
total thermal conductivity 
frozen thermal conductivity 
molecular weight of species i 

shape parameter, n = 1 for axisymmetric and n = 0 for two-dimensional 
total number of moles 
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pressure 

total Prandtl number 

heat -transfer rate to the wall 

y component of the heat -flux vector 

body coordinate, shown in sketch, appendix A 

temperature 

velocity in direction of boundary -layer flow (i.e., velocity in x 
direction) 

free -stream velocity 

velocity normal to direction of boundary -layer flow (i.e., velocity in 
y direction) 

coordinate along wall 

mole fraction of species i 

coordinate normal to wall 

pres sure -gradient parameter defined in equation (A21) 
transformed coordinate defined by equation (A12b) 
viscosity 

transformed coordinate defined by equation (A12a) 

mass density 
PM- 

PwMw 


Subscripts 

value at onset of dissociation or ionization 

boundary -layer -edge value 

stagnation -point value 

reference value, see table I 

wall value 

free -stream value 
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PRESENTATION OF EQUATIONS 


The problem considered is that of convective heat transfer in various 
gases. The equations necessary for describing this phenomenon are the conser- 
vation of mass, momentum, and energy equations. Along with these, relations 
describing the fluid state and transport properties are required. These equa- 
tions are developed in appendix A and the resulting equations are presented 
below. 


Heat -Flux Equation 

The wall convective heat -flux equation for a chemically reacting mixture 
of gases in thermochemical equilibrium is 


k w /dhN _ Mw ( dh\ 
c Pw W/w Pr w W/ w 


( 1 ) 


This equation can be solved when the transport properties and enthalpy gradi- 
ent at the wall are known. For this study, the unknown in equation (l) is the 
enthalpy gradient which must be obtained by solving the boundary -layer 
equations . 


Similar Boundary -Layer Equations 

An appropriate set of equations for studying the effect of gas composi- 
tion on the wall enthalpy gradient is the similar form of the boundary-layer 
equations. These particular equations are exact when similarity holds, such 
as in the stagnation region of a body, and they allow approximate prediction 
of the heating -rate distribution over bodies. 


The similar boundary -layer equations in the familiar | and rj coordinate 
system (see eq. (A12) ) may be written in the following form: 


(cpf")’ 


+ ff" + 3 



0 


( 2 ) 


where 



+ fg* 





0 
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( 3 ) 
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ar}.d where the prime superscript represents differentiation with respect to rj . 
The Prandtl number in equation (3) is the total Prandtl number obtained when 
the total specific heat and total thermal conductivity are used. 

The boundary conditions on these equations are as follows: 

f(o) = 0 f'(oo) - 1 
f * (o) = 0 

g(o) = g w g(“) - 1 

The thermodynamic and transport property terms in these equations are 
treated as known functions of the static enthalpy. The pres sure -gradient 
parameter (3, wall temperature, pressure, and the term Ue 2 /H s are treated 
as parameters. 


Equations for Heating Rate and Heating-Rate Distribution 

Solutions to equations (2) and (3) are in the g and tj coordinate system; 
therefore, with the aid of the transforming equations in appendix A, equa- 
tion (l) is rewritten as 


Hw “ 


r Pw^w u e^s 

Pr w JJ[ 


g* ( 0 ) 


(4) 


To obtain the heating-rate distribution over a body it is convenient to 
normalize equation (4) by the heating rate at a stagnation point given in 
reference 2 as 


1w 0 




M-w 0 



(5) 


Using equations (4) and ( 5 ) and assuming an isothermal surface, at a tempera- 
ture small compared to the stagnation temperature, we may write the following 
heating-rate distribution equation 


1 

II 

1 

n 

r 

( Pw 4 / u e 4 

V W 0 J vv g* (0) 

q w 0 

J 


/u e \ sn 

{d r 

“dCue/U^)" 

dx 

sh°) 

0 


( 6 ) 
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THERMODYNAMIC AND TRANSPORT PROPERTIES 


Solutions to the boundary -layer conservation equations depend on a 
knowledge of the thermodynamic and transport properties of the gas in ques- 
tion. The thermodynamic properties of pure gases can be calculated to good 
accuracy from spectroscopically determined constants. The situation concern- 
ing transport properties is more uncertain and their calculation depends on 
assumptions concerning the intermolecular potentials at the lower temperatures 
and on many uncertainties regarding collisions and diffusional phenomena at 
the higher temperatures. Despite these uncertainties, properties so calcu- 
lated are probably representative of the actual values and, at least, should 
show how differences in gas composition affect heat transfer. This study was 
confined to the gases, air, nitrogen, hydrogen, carbon dioxide, and argon. 

The actual values of the properties for the various gases were compiled from 
information in references 6 to 15 . 

The properties were curve fitted as a function of enthalpy and used in 
the numerical solutions to the boundary -layer equations. This is explained in 
appendix B and the coefficients of each curve fit are tabulated in tables II 
through V. 

To use these properties the gas considered must be in local thermochemi- 
cal equilibrium and its atomic composition must correspond to that given by 
assuming its initial state was at standard conditions of pressure and tempera- 
ture. Hence, using these properties in the boundary -layer equations implies 
equal diffusivity among the various species at all points in the boundary layer. 

Before discussing the final results of the study, it is of interest to 
examine the variation of gas properties with enthalpy as these quantities 
enter the basic conservation equations as coefficients. 

The gas-density ratios p e /p are plotted against the enthalpy ratio 
h/H s in figure 1. This figure is representative of the change in density 
with enthalpy at the stagnation region of a blunt body for a flight velocity 
of 30,000 ft/sec and a pressure of 0.1 atm. A significant difference in den- 
sity between the various gases occurs at intermediate values of the enthalpy 
ratio. However, as stipulated in reference 1, differences of this magnitude 
probably have little influence on the wall enthalpy gradient obtained from the 
solution to the boundary-layer equations. Further investigation was under- 
taken to verify this point and the results are discussed later. 

Figure 2 presents the variation of pfi/p w |U w (hereafter defined as cp) 
with h/H s for each gas. Air, nitrogen, carbon dioxide, and argon behave in 
a similar manner, showing some differences in level with particular values of 
h/Hs* Hydrogen exhibits the smallest variation of cp. ■■ 

Figure 3 shows the changes in cp/Pr with h/Hs for each gas. This term 
reflects the changes in the ratio k/cp superimposed on the density variation 
seen in figure 1. Fluctuation of the cp/Pr curves from h/H s = 1 to 
h/H s = 0 is related to the chemical reactions that take place during the 
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Figure 3 .- Variation of cp/Pr with enthalpy for various gases. 


recombination of various species through the boundary layer. These quantities 
depend both on the concentration gradients and concentrations of each species 
and therefore assume rather complicated behaviors. The fluctuations of cp/Pr 
are usually associated with the completion and onset of various reactions. 

For example, the species concentration of C0 2 taken from reference 12 and 
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plotted in figure 4 shows that the first fluctuation occurs during the 
formation of the maximum amount of CO from free 0 and C and that the next nota- 
ble fluctuation occurs during the final formation of C0 2 from the various dis- 
sociated species. The effect of this cp/Pr variation on the solution to the 
boundary-layer equations will be discussed later. 



h/H s 


Figure 4.- Gas composition of C0 2 (mole per original cold mole of C0 2 ) . Variation with static enthalpy 

(taken from ref. 12). 
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Pressure level can affect the property variations of each gas as 
illustrated in figure 5 where cp/Pr for N 2 and C0 2 for two pressures is s 
For N 2 there are only small differences in cp/Pr for the two pressures, 
whereas for C0 2 there is a marked difference. The other gases do not exhi 
such wide differences with pressure level as C0 2 , but the need to investig 
the effect of pressure changes on the boundary-layer solutions for all the 
gases is evident. 




• The comparisons made above are inconclusive insofar as telling how the 
differences in gas properties might affect solutions to the boundary -layer 
equations and, in particular, the wall enthalpy gradient. The next step is to 
solve the differential equations for various values of pressure, pressure gra- 
dient, boundary -layer -edge velocity, and wall temperature. The range for 
these was generally chosen as follows: 

0 < 0 < 1 

500° R < T w < 1850° R 


10 " 3 < p e < 10 atm 
RESULTS AND DISCUSSION 


Similar Solutions 

Solutions to the boundary -layer equations for each of the gases when 
certain external flow parameters and wall conditions were varied are summa- 
rized next. Generally, each gas behaved in the same way; for example, the 
heat -transfer parameter g*(o)/l-g w increased with increasing wall tempera- 
ture and increasing pressure -gradient parameter 3; it decreased with increas- 
ing flight velocity or total enthalpy; for practical purposes it was unaffected 
by changes in the dissipation parameter u e 2 /Hs and the pressure level. 

A correlation which included all the input parameters was needed to illus- 
trate the important results of the solutions. Such a correlation was pre- 
sented by Pay and Riddell (ref. 2) and subsequently by Kemp, Rose, and Detra 
(ref. 3). These investigators found that for dissociated air (for flight 
velocities to 30,000 ft/sec) the heat-transfer parameter g* (o)/l-g w corre- 
lated with cp e . This term reflects implicitly the changes in flight velocity 
and wall temperature; that is, at constant wall temperature it decreases with 
increasing flight velocity; at constant flight velocity it increases with 
increasing wall enthalpy or temperature. In addition cp e appeared to account 
for the over-all changes in thermodynamic and transport properties through the 
boundary layer, implying that subsequent differences in the numerical values 
of the transport properties would not affect the correlation equation. The 
investigators of reference 3 found that the pres sure -gradient parameter was 
accounted for by assuming that the enthalpy gradient varied as (l + constant N /p) 
and also that the dissipation term (u e 2 /Hs)[(cp - cp/Pr)f , f"]' affected the 
solutions to a minor degree. This correlation was attempted for each of the 
gases considered herein. 
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The results are presented in figures 6 through 10. In these figures the • 
heat -transfer parameter g* (o)/l-g w is divided by Pr w and by a term vhich 
accounted for the pres sure- gradient parameter (i.e., (l + constant s/T) ) • The 
various symbols represent different values of the pressure -gradient parameter. 


(l-g w ) Pr w 


= .667(l+.096v^)(^ e ) 0 475 
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Figure 6.- Heat -transfer parameter correlation for air. 
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Figure 7 • " Heat -transfer parameter correlation for nitrogen. 












Figure 8.- Heat-transfer parameter correlation for hydrogen. 


Symbols are flagged to represent solutions for different wall temperatures and 
different pressures. Symbols are filled to represent solutions for different 
values of u e 2 /H s . For the investigated range of parameters the solutions 
correlate about straight lines of the same slope fitted to the solutions for 
constant values of u e 2 /H s . The two straight lines shown on each figure 
represent curve fits through the solutions for u e 2 /H s = 0 and u e 2 /Hs = 0.5« 

The first important aspect of the correlation is the relatively small 
dependence of the heat-transfer parameter on f3« The term (l + constant ./p) 
found in the ordinate of each figure is never much larger than 1.0 because 
the value cf the constant, although slightly different for each gas, is numer- 
ically small. 

The correlation figures point out that the dissipation term also has a 
relatively small effect on the enthalpy gradient. This is clearly demon- 
strated in figure 9 where the results for C0 2 are plotted. At three distinct 
locations along the correlating line for u e 2 /H s = 0 the dissipation parameter 
was varied between 0 and 1 for a constant pressure and wall temperature and a 
solid line (practically horizontal) joins these solutions. For each succes- 
sive change in u e 2 /H s the value of T e increases because the static 
enthalpy decreases at the boundary-layer edge while the enthalpy gradient 
changes a small amount. The largest changes in the enthalpy gradient occur at 
the largest values of T e > but for practical purposes these are negligible. 
Each of the gases showed a similar dependence on the dissipation parameter. 
This allows a convenient equation for the enthalpy gradient to be expressed in 
terms of the stagnation value of T e and (3. These equations are presented in 
each of the figures. The advantage of these equations is that the enthalpy 
gradient ratio in the heating -rate distribution (see eq. (6)) can be deter- 
mined without computing the local value of the density-viscosity product. 
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Next, these correlations demonstrate the effect of changes in vail 
temperature and pressure. As the vail temperature increases, cp e increases. 

At the same time, the heat -transfer parameter increases so that a straight 
line joins solutions for identical values of u e 2 /H s . Note that differences 
in Pr w , resulting from changes in vail temperature, did not introduce scatter 
into the correlations. Likevise, changes in pressure cause changes in cp e 
vith corresponding changes in the heat -transfer parameter such that the cor- 
relation line is maintained. The scatter of the results about the correlating 
line due to pressure and temperature changes is small for practical applica- 
tions. These changes introduce more scatter for C0 2 (fig. 9) than for air, 
hydrogen, and nitrogen. An estimate of the magnitude of this scatter is 
obtained by comparing some of the results for u e 2 /H s - 0 vith the straight- 
line curve fit to these points. It can be seen that changes in pressure level 
introduce more scatter than vail temperature. 



Figure 9*" Heat-transfer parameter correlation for carbon dioxide. 


Argon (fig. 10) did not correlate as well as the other gases. For this 
gas, only the solutions for the lowest value of the wall temperature were used 
for obtaining the straight-line fits. At each velocity (constant peMe) 
increasing the wall temperature increases both the enthalpy gradient and cp e . 
However, decreasing cp e by maintaining a constant wall temperature and 
increasing the flight velocity does not cause a corresponding decrease in the 
enthalpy gradient. This is in direct contrast to the solutions obtained for 
the other gases. As a result, the slope of the correlating line is much 
smaller than for the other gases. Pressure also affects the solutions for 
argon more than other gases ; for example, at 20,000 ft/sec, results for 10' 3 
and 10 ~ 1 atm agree very well but differ from that for 10 atm by about 10 per- 
cent, although this is not considered significant for practical applications. 


14 




+. 


Figure 10.- Heat -transfer parameter correlation for argon. 


The behavior of argon is quite similar to that of air when the total enthalpy 
reaches the value where ionization significantly affects the viscosity used in 
evaluating Te (see, e.g., ref. 16 ) . 

Correlations in terms of Te are not without shortcomings. First, the 
correlation equations can he used only over the range of Te for which they 
are developed. This is substantiated in the results for air given in refer- 
ences 1 6 and 17 and will be demonstrated later in this study. Secondly, a 
single correlation equation in terms of Te is not obtained for all the gases 
and therefore the results are not general enough for computing the heat trans- 
fer in gases other than those being considered. Finally, the correlation 
lacks the quality of accounting for the constancy of the enthalpy gradient 
with varying u e 2 /Hs. Two unsuccessful attempts were made to eliminate these 
shortcomings. The integrals of cp and cp/Pr across the boundary layer were 
investigated as single correlating factors for all gases. However, they 
resulted in individual correlations for each gas which were no better than 
those for which Te alone was used. Hence, a detailed study of the effect of 
the property variations on the solutions was undertaken. The results are 
presented next and are then used to obtain a single correlation which is use- 
ful for estimating heating rates in various gases. 
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Effect of Gas Properties on Similar Solutions 


Some of the solutions in carbon dioxide gas are considered first since 
they are representative of the results for other gases. Variations of veloc- 
ity f ! , enthalpy g, and the property terms across the boundary layer are 
given in figures 11, 12, and 13 for a pres sure -gradient parameter of l/2 and 
for three flight velocities. Two values of the dissipation parameter were 
used, one representing a stagnation point and the second, a point on the vehi- 
cle where the local boundary-layer outer-edge velocity is ^ H s /2. The local 
boundary-layer -edge velocity influences the solutions in two ways: first, it 

changes the static enthalpy distribution across the boundary layer which in 



0 1234 0 1234 



0 12 3 4 

V 

Figure 11.- Comparison of boundary -layer profiles; C0 2 , =■ 10,000 ft/sec, T w - 500° R, and 

p - 0.1 atm. 
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turn changes the fluid property variations; secondly, it introduces the 
dissipation term in the energy equation (see eq. (3))- It will he shown sub- 
sequently that the latter effect is very small and therefore the comparison of 
the pairs of solutions in figures 11 through 13 shows the effects of varying 
property distributions through the boundary layer. 


In figure 11, the enthalpy, velocity, and property terms vary smoothly 
between their wall and boundary -layer -edge values . The relatively small 
change in property profiles, brought about by including u e 2 /H s , does not 
affect the enthalpy gradient at the wall to any significant extent. In fig- 
ure 12 where the total enthalpy has been increased, significantly different 



0 1234 0 1234 

T) r) 



Figure 12.- Comparison of boundary -layer profiles; C0 2 , = 20,000 ft/sec, T v 

p = 0.1 atm. 
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property variations are obtained, and yet the enthalpy profiles, especially 
near the wall, do not reflect these differences. This is especially true for 
the enthalpy gradients at the wall which differ for this case by less than 
10 percent. Similar conclusions are reached from figure 13. The results 
indicate that property variations "far away" from the wall (especially for 
cp/Pr) have a rather small effect on the enthalpy gradient at the wall and 
therefore on the heat transfer. Another interesting aspect of the comparison 
which will be used later is the behavior with total enthalpy (or flight veloc- 
ity) of the cp and cp/Pr terms near the wall. The absolute value of the slope 
of these terms near the wall increases with flight velocity while the wall 
enthalpy gradient decreases (recall the results for varying free -stream veloc- 
ity in figs. 6 to 10). 



v 



0 12 3 4 
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Figure I 3 .- Comparison of "boundary -layer profiles; C0 2 , 

p = 0.1 atm. 


0= 5, u e 2 /H s =0 

13 - .5, u e 2 /H s = .5 


Uco = 30,000 ft/ sec, T w 


500° R, and 
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♦ It was mentioned above that the effect of including the dissipation term 
in the energy equation became evident only through its effect on the fluid 
property variation through the boundary layer. To show this, the magnitude of 
various terms in the energy equation (see eq. (A26) ) for the flight velocity 
of 20,000 ft/sec is compared in figure 14 with those obtained when the 



v 

Figure 14. - Magnitude of terms occurring in energy equation; C0 2 , U OT - 20,000 ft/sec, T w = 500° R, and 

p = 0.1 atm. 

dissipation term was zero. The value of U e 2 /H s was one -half in both cases. 
This figure shows that the dissipation term itself (u e 2 /H s ) [(cp - cp/Pr)f *f "] * 
has little influence on the solution and especially on the wall enthalpy gra- 
dient. Spot -check solutions for other velocities and for other gases gave the 
same result. 

It is also appropriate at this point to indicate how the density term 
p e /p affects the determination of the enthalpy gradient. Several solutions 
to the boundary-layer equations for different density distributions showed 
that the solutions were insensitive to arbitrary variations of p e /p. For 
example, using the p e /p distribution for C0 2 and the cp and cp/Pr distribu- 
tion for N 2 to obtain a stagnation -point solution at 30,000 ft/sec resulted in 
a negligible change in enthalpy gradient from the corresponding nitrogen 
solution. 
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Thus far, it has been shown that the variation of the terms cp and cp/Pr , 
"close" to the wall rather than that "far" from the wall affects the enthalpy 
gradient and that the rate of change of these property terms with rj near the 
wall influences this gradient. For this reason a single correlation curve 
with cp e is not obtained for all gases. This is illustrated in figure 15 
where the stagnation-point profiles for some of the gases are plotted for a 
flight velocity of 30,000 ft/sec, a pressure of 0.1 atm, and a wall tempera- 
ture of 500° R. Also included are the computed values for g' (o)/l-g w . 

Observe that the values of cp e for the gases differ but that the enthalpy 
gradients do not, provided the variations of the property terms near the wall 
are similar (i.e., compare the enthalpy gradients and the property profiles 
"near" the wall for C0 2 , N 2 > and A) . These results suggest that any general 
correlation should probably be based on properties evaluated "near" the wall. 





V 



0 12 3 4 

V 


Figure 15 • - Boundary -layer profiles for various gases at a stagnation point; - 30,000 ft/sec, 

p 0.1 atm, and T v - 500° R. 
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• The insensitiveness of the solutions to the variation of the transport 
properties in the higher energy portions of the boundary layer indicates that 
precise evaluation of the transport properties at the higher enthalpies may 
not be critical for determining the equilibrium convective heat transfer. 
Calculations were made to determine how much the properties could be changed 
before significant differences in the enthalpy gradients appeared. Typical 
results are presented in figure 1 6 where the profiles for a stagnation point 
in carbon dioxide are plotted. First, the true cp/Pr variation was changed 
by fixing the Prandtl number constant (at its wall value) across the boundary 
layer. The velocity and enthalpy profiles in the lower energy portion of the 
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boundary layer are not changed and the enthalpy gradients differed by less 
than 5 percent. Next, solutions for several different cp profiles were 
obtained by arbitrarily fixing cp at a constant value at some point in the 
boundary layer; the enthalpy and velocity profiles and the enthalpy gradient 
did not change significantly. (For these cases the cp/Pr variation with 
constant Prandtl number described above was used.) Calculations for various 
velocities demonstrated that large differences in p e or (k/cp) e do not sig- 
nificantly affect the value of the enthalpy gradient. One additional case 
was computed that is worthy of comment. For a stagnation -point solution at 
30,000 ft/sec, cp/Pr (where Pr = Pr w ) was allowed to vary normally until t) 
reached a value of 1; from there cp/Pr was increased until at the boundary- 
layer edge it was 25 times its normal value. The enthalpy gradient was 
increased only 50 percent over the correct value. Concurrent with the present 
investigation, reference 18 determined the effect of uncertainties in the 
thermal conductivity of air on the stagnation -point heat transfer. Their con- 
clusions are in substantial agreement with those made above. These results 
vividly illustrate that correlation equations in terms of boundary -layer 
outer -edge properties should be used only over the range of cp e for which 
they are derived. 


Correlation of the Heat -Transfer Parameter in Terms of 
Low Temperature Properties 


The results described above were used as a basis for obtaining a single 
correlation of the wall enthalpy gradients for all the gases. It seemed rea- 
sonable that the wall gradient would correlate on an "average" derivative of 
the "near-wall" property terms, and during the course of making arbitrary 
changes in the property variations it was determined that changes beyond the 
point where the cp/Pr variation reached a minimum or inflection did not seri- 
ously alter the enthalpy gradient at the wall. The correlation was obtained 
by using the absolute value of the average slope of cp with h/H s , that is, 


dcp 

d(h/H s ) 



dcp 

_d(h/H s ) _ 

hp ^w 
H s " H s 




Both hp and cpp were evaluated at the point in the boundary layer where dis- 
sociation (or ionization for argon) was just beginning. In particular, this 
point was chosen so that for the dissociating gases Z (the ratio of molecular 
weights) was 1.01 and for argon a (the degree of ionization) was 0.001. The 
points in the boundary layer where these conditions prevailed were associated 
with the points where cp/Pr obtained its first minimum or inflection (e.g., 
see cp/Pr in fig. 15) . At the enthalpy hp, only small amounts of dissoci- 
ated or ionized species are present in the gases, but the concentrations and 
concentration gradients of the various species are sufficient to cause the 
ratio k/cp to remain constant or decrease as reflected by the cp/Pr term. 

Figure 17 shows the correlation of [g'(o)/l-g w ] Jp w p w /(pp) T= ^ 00 o R - The 

stagnation-point solutions of the present investigation are included (with the 
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H s 

Figure 17 . - Correlation of stagnation -point solutions for all gases. 


exception of some hydrogen cases for which no dissociation occurred) and are 
found lo correlate within ±1 6 percent. The term j Pw^wAppOij^oo 0 R was 

included to reduce some scatter due to wall temperature variations. Also cor- 
related are several solutions of reference 2 (Lewis number equal to 1.0) with 
values of hp and cpp given in reference 6 corresponding to a pressure of 
0.1 atm and Z (ratio of molecular weights) of 1.01. The solution described 
earlier where (k/cp) e was 25 times the true value is also included. For this 
case hp and cpp were evaluated at the point in the boundary layer where cp/Pr 
reached a minimum. The equation describing the correlation is 


g 1 (o) / Pw^w 

1 " g w J (^1=500° R 


0.324 


1 ~ 7p 

L(h D /H s ) - (h w /H s )J 


-O . 15 


(7) 


If the enthalpy gradients for other values of p and u e 2 /H s are desired, they 
can be obtained by multiplying the corresponding stagnation-point gradient by 
the factor 0.93 (l + 0.10J13). This factor uses an average value of the con- 
stants in the .expression (l + constant JF) in figures 6 to 10. 


Other choices for determining the slope of the "near -wall" properties 
were also investigated. The best results were obtained by using a weighted 
average slope of cp/Pr with h/H s over the same limits of integration used 
above. This resulted in a correlation of the stagnation-point solutions to 
within ±13 percent. This is about the smallest percentile range that can be 
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obtained using property derivatives with respect to h/H s as a correlation * 
parameter. The equation describing this result is 

-O . 12 


(8) 


Obviously, the simplicity in applying equation ( 7 ) outweighs the small corre- 
lation improvement given by equation (8). 

The correlation in terms of the property derivatives does not contain the 
shortcomings inherent in the correlations in terms of T e • First, a single 
correlation equation arises which can be used to predict the enthalpy gradient 
in gases other than those being considered. In application, this becomes a 
rather simple task since the correlation depends on relatively low temperature 
properties which are calculable for pure gases or gas mixtures by the methods 
described in reference 19* The correlation can be extrapolated into a range 
of enthalpy where both dissociation and ionization occur together in the 
boundary layer, as will be shown subsequently by demonstrating reasonably good 
agreement between heating rates calculated from the present correlation equa- 
tion and those predicted by theories accounting for the effects of ionization 
on the transport properties and with actual shock -tube data. Finally, the 
correlation accounts for the constancy of the enthalpy gradient with u e 2 /ll s . 


Pw^w 


g*(p) r 

^ ~ J (pOt^oo 0 r 


= 0-37 


(9/Pr)^ - (cp/Pr)J 


'K/Ks 


(9/Pr) d(h/H s ) 


Stagnation -Point Heat Transfer 

The heat-transfer parameter, g'(o)/l-g w , was used to describe the results 
of the similar solutions to the boundary -layer equations . The next concern 
is how these results apply to the determination of the heat transfer to a body. 

First, it is informative to develop specific equations for stagnation- 
point heat transfer in each gas so that the relative levels of heating rate 
for comparable conditions of flight velocity and stagnation -point pressure can 
be observed. 

If a Newtonian pressure distribution in the vicinity of the stagnation 
region is assumed, the equation for the velocity gradient in equation (5) 
becomes 



(9) 


For simplicity p^ is neglected in what follows. Substituting equation (9) 
i nto e quation (5) and multiplying both sides of the resulting equation by 
n/R/Po gives 
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H s /Q ( PyH-Hf) 2 

Pr w p 0 Pe Q 


g’ (o) 


( 10 ) 


The right side of equation (10) was computed for each gas for a range of 
stagnation pressures from 10“ 3 to 10 atm, wall temperatures from 500 ° to 
1 , 850 ° R, and total enthalpies corresponding to flight velocities from 
10,000 to 30,000 ft/sec. In this determination the value of g* (o) from the 
■boundary -layer solutions was used and not the correlations described above. 
The results were correlated within 10 percent by an equation of the form 




(ID 


where is dimensionless and equals U^, divided by 10,000 ft/sec. In most 

cases the deviation of the computed value from that given by equation (ll) was 
much smaller than 10 percent. Equation (ll) is plotted in figure 18 and the 
constants G and N are tabulated. For a given flight velocity, argon has the 
highest stagnation -point heat transfer followed successively by C0 2 , air, N 2 , 
and H 2 . For a given flight irelocity, the stagnation -point heating in C0 2 , 
air, and nitrogen is nearly the same. For hydrogen, the heating rate is not 
very large compared to the other gases over the range of velocity considered. 

Although equations of the form given by equation (ll) are very useful, 
their extension to higher velocities and to other gases is not straightforward 
since the constants C and N are valid only in the range of velocity consid- 
ered here and are not readily expressed as functions of gas properties. 
Therefore it becomes necessary to develop a more general equation. Substi- 
tuting equation ( 7 ) into equation ( 5 ) results in the following heating -rate 
equation for an axisyrametric body 


% 


= 2^24 I, 


Pr w n / 


2 (pM-) t= 500° r 


clu ( 

dx 


1 - <PD 


-,-0.15 


L(h D /H s ) - (h w /H s )J 


Hsd-Sw) (12) 


Furthermore, using equation ( 9 ) and the approximation, H s » U, C0 2 /2, we can 
rewrite equation (l2) as 


_ 0,360 1 [( p|i ^T=500 o r1 ( 1 - ° ' 15 

1 - Sw = Pr w 4 J 2 P O Pe 0 \ h D ' h J 

l/2 

To obtain equation ( 13 ) in the units, Btu/ft 2 ~sec (ft/atrn) , 
rewritten 


1.7 


it is 


( 13 ) 
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Figure 18.- St agnation -point heat transfer for various gases. 
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sec, and 


sed in slugs/ft 3 , p in slugs/ft sec, h in ft 2 / sec 2 , 
in atm. Equation (l4) can be used to estimate the heating 


rate in various gases. To apply the equation, knowledge of the low tempera- 
ture properties of the gas along with the stagnation -point pressure and den- 
sity is required. 

To substantiate the theoretical predictions made above, a comparison of 
theory and experiment is given in figures 19(a) to (d). Here, heating-rate 
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Figure 19*- Continued. 


data from various sources are plotted against (H s - h w ) . For air, the shaded 
area represents the majority of shock-tube data from a large number of sources. 
The solid lines represent equation (ll) with the appropriate values of C and 
N. The data and equation (ll) agree well for all the gases. Also, equa- 
tion (ll) agrees very well with the theory of reference 4- for air and CO2 and 
reference 20 for air. 

Equation (l4) is also shown in figure 19 evaluated for p 0 = 0.1 atm. 

This equation, as expected, agrees well with equation (ll) up to 30,000 ft/sec. 
To demonstrate the utility of this equation, which uses the correlation shown 
in figure 17, it was extrapolated to substantially higher enthalpies. The 
agreement with the data is good. Also, the theory of reference 20, which 
includes the use of the high temperature transport properties of air, and the 
extrapolation of equation (l4) agree favorably in the extrapolated region 
where ionization and dissociation occur together in the boundary layer. There- 
fore, it appears that the heating rate in various gases can be predicted ade- 
quately to very high enthalpies by an equation based on relatively 
low-temperature properties. 

Next, the problem of heat distribution on the surface of the body is 
considered. 
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Figure 19- - Concluded. 


Heating-Rate Distribution 

Equation (6) shows that the distribution of heat transfer around a body 
depends explicitly on the inviscid flow as reflected in the boundary-layer- 
edge velocity and surface pressure and implicitly on the inviscid and viscid 
flow as reflected in the variation of g' (o)/g^(o) with pressure gradient. An 
investigation of the magnitude of this ratio of enthalpy gradients is pre- 
sented next and it shows that including this ratio in the heating -rate distri- 
bution equation modifies the results by a small amount. 

Figure 20 presents the variation of g' (o)/g£(o) for an axisymmetric body 
with pres sure -gradient parameter for each of the gases for a 'single value of 
the dissipation term u e 2 /Hs = l/2 at three flight velocities and at a speci- 
fied wall temperature and pressure. Each gas shows an increase in this ratio 
with increasing pressure gradient. Positive values of (3 imply favorable 
pressure gradient. The ratio of enthalpy gradients for any particular gas 
varies by about 10 percent over the range of (3 considered. It is noteworthy 
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Figure 20,- Enthalpy-gradient variation with pres sure -gradient parameter for ax i symmetric configurations. 


that the curves do not always pass through 1.0 at {$ = l/2 as a direct conse- 
quence of including the dissipation term in the equations and, as pointed out 
hefore, its effect on the solution is small. For all values of the dissipa- 
tion term used in this investigation, the effect never exceeded 10 percent. 

If the dissipation term is neglected, the correlation equations presented in 
figures 6 through 10 can be used to express the axisymmetric value of 

g'(°)/go(°) as 

g'(o) i + c./ir ( 15 

g’(o) " 1 + 0.707c 

where c is the constant given in figures 6 to 10. 
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- It is interesting to compare the variation of g , (o)/g ( ’(o) with that 
obtained in reference 21 for an ideal-gas solution to the boundary -layer equa- 
tions with cp and Pr set equal to 1. The dashed curve in figure 20 represents 
such an ideal-gas solution, and as shown, this approximation represents quite 
adequately the real-gas variation. This ideal solution is also compared with 
equation ( 15 ) for carbon dioxide (see = 10,000 ft/sec) and the comparison 
is very good. Either method for obtaining g* (o)/g£(o) is adequate for the 
range of pressure -gradient parameters considered. Ideal-gas solutions to 
pres s\ire -gradient parameters of 1+.0 are tabulated in reference 21. 

It is concluded that the transport and thermodynamic property variations 
play a minor role in determining the heating-rate distribution, at least for a 
practical range of the pres sure -gradient parameter, whereas the external flow 
field is of major importance. 


CONCLUSIONS 

An investigation of the effects of gas composition on the equilibrium 
convective heating rate and heating -rate distribution resulted in the follow- 
ing conclusions: 

1. The effect of the transport property variations in the higher energy 
portions of the boundary layer was not reflected in the wall enthalpy gradient. 
At large ij changes by a factor of 3 in thermal conductivity or viscosity 
changed the wall enthalpy gradient and therefore convective heating rate less 
than 10 percent. For the range of parameters investigated, this indicates 
that precise evaluation of transport properties at high temperatures is not 
critical to equilibrium convective heat -transfer calculations. 

2. The wall enthalpy gradient obtained f jTOIii SOlut ions to the similar 
form of the equilibrium boundary-layer equations for air, nitrogen, hydrogen, 
carbon dioxide, and argon can be correlated on low temperature properties. 

This correlation can be used to extend the present results to other gases and 
to higher total enthalpies. 

3- The stagnation -point heat transfer depends on gas composition. For 
the same body radius, total pressure, and flight velocity, argon gives the 
highest heat transfer and hydrogen gives the lowest value. Air, nitrogen, and 
C0 2 give about the same intermediate value of heat transfer. 

4. The heating-rate distribution on a body is affected to a minor degree 
by the gas composition. For practical application, the heating -rate distribu- 
tion to blunt bodies can be obtained from inviscid flow considerations alone. 
Refinements to this can be obtained by using solutions to the low speed form 
of the boundary-layer equations. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Feb. 23, 1965 
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APPENDIX A 


DERIVATION AND SOLUTION OF EQUATIONS 


y 



The analysis is restricted to a gas in thermo chemical equilibrium whose 
chemically reacting species are considered a mixture of perfect gases. For 
such a gas, the general equations of change for a body-oriented coordinate 
system (see sketch above) subject to Prandtl's boundary -layer assumptions are: 


Bx 


pu |s. + p T |u + a 

ox ay dx 
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by by 
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(- I )- 0 

(A2) 

Su 2 /2 n 
by / 


(A3) 


where n = 1 for axisymmetric bodies and n = 0 for two-dimensional bodies. 
Equation (Al) is the continuity equation and represents the sum of the indi- 
vidual species equations. Equation (A2) is the momentum equation which 
relates the momentum change to the pressure and viscous shear forces acting on 
the system. Equation (A3) is the energy equation which relates the enthalpy 
change to heat addition by viscous stresses and the normal component of the 
heat -flux vector, qy. 

The normal component of the heat -flux vector (neglecting thermal diffu- 
sion and radiation) is related to the temperature and molar concentration of 
the gas by the equation (see ref. 19 ) 
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The first and second terms on the right side of equation (A4) represent, 
respectively, heat addition due to conduction and due to diffusion of species 
across the boundary layer. Enthalpy is introduced into equation (A4) as fol- 
lows. The molar concentration and enthalpy of the gas in equilibrium vary 
with temperature only (the pressure is considered constant across the boundary 
layer) so 


and 
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and by definition the total specific heat is 
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= X Cp i Gi 
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Therefore (A4) , with the aid of equations (A5) and (A 6) and the above defini- 
tion of specific heat, becomes 
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where the bracketed term in equation (A8a) is the so-called "total" thermal 
conductivity k. 

Now, the energy conservation equation (eq. (A3)) is rewritten with the 
aid of (A8b) as 
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where Pr represents the total Prandtl number of the gas obtained when total* 
values of thermal conductivity and specific heat are used. Equation (A9) may 
be rearranged to give the energy equation 


cu & + pv an . I- (-as + , 0 

Bx Sy dy |_ \dy Pr dy/_ 


Equations (Al), (A2), and (A10) are the conservation equations for a reacting 
gas in thermochemical equilibrium. The solution to these equations is 
required when the heat flux to. a wall is computed. This is apparent by 
inspection of equation (a 8) (evaluated at the wall), 


dh\ _ Mw f dh\ 


W C P W w/ w Pr w w/ w 


To put equations (Al), (A2), and (A10) in a form more suitable" for solu- 
tion, the following transformations are introduced: 


I = / PwPw u e r dx 


The transformed differential operators then become: 


PwPw u e r 
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It is convenient to make the following definitions : 
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The continuity equation is satisfied by the introduction of the stream function 
t defined by equations (Al8a) and (Al8b) . Applying the above definitions and 
the transformed differential operators results in the following form of the 
equations : 
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(A19) 


2|(f n g| - f|g n ) (A20) 


where the subscripts denote partial differentiation with respect to that 
variable and 


3 = 2 


d In u e 
d In | 


(A21) 


Equations (A19) and (A20) comprise a set of simultaneous, nonlinear partial 
differential equations. The local similarity assumption is used to solve this 
system of equations. In this approach, the terms on the right side of equa- 
tions (A19) and (A20) are assumed small compared to other terms in the equa- 
tions so that methods applying to solutions of ordinary differential equations 
can be applied. This means that changes of the dependent variable with i 
are negligible and that the quantities at the boundary-layer edge assume their 
local values as specified by external flow conditions. Further discussion of 
local similarity and its application can be found in references 1 and 3* 


The similar form of transformed boundary -layer equations for the conser- 
vation of momentum and energy form a coupled pair of ordinary, nonlinear dif- 
ferential equations of third order in f and second order in g. These are 
written as 
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where the prime denotes differentiation with respect to r), that is, d/dq. 

These equations comprise an initial value problem requiring five conditions 
at T) = 0 (namely f(o), f T (o), f n (o), g(o), and g T (o)) for a unique solu- 
tion. Since only three of these conditions are considered known, one must 
seek a solution using iterative techniques. This is accomplished by imposing two 
additional conditions on the equations as r\ °°. The boundary conditions for 
the equations in this case are: 

f(o) =0 f («>)-* l 
f' (o) - 0 

g(°) = g w g(°°) - 1 

It is advisable to eliminate the explicit forms of the property deriva- 
tives from the equations since the accuracy of these terms would be highly 
questionable at high temperatures. This can be done by introducing the fol- 
lowing integrating factor 
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Equation (A22) may be rewritten as 
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which when integrated from 0 to q and solved for f ,r gives: 
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Equation (A23) when integrated from 0 to rj and solved for g* /Pi’ gives 
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Then the conservation equations may he written as : 
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Equations (A27) through (A3l) along with 


f = 



(A32) 


are the six equations to he solved simultaneously to determine f^, g^, and 
profile distributions throughout the boundary layer. 


The iterative technique that is used to obtain a solution is one of 
successive approximations based on the Newton -Raphson method (see ref. 22). 
The functions, f T and g, are expanded in a Taylor series where q -» oo and 
the higher order terms neglected. Thus, for f*(°°) = F(f^), g^) and 

g(°°) = Sw) we have 
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By approximating the differentials with finite differences we may obtain 
improved approximations for f” and gA. as follows: 
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where the subscript o represents the corrected values so that f^(°°) = 1 and 
g»(oo) - i } and we have: 
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(A 38 ) 


The iteration scheme is started with two initial guesses for f^. and g' 
(denoted by f-w^j fw 2 > gw±> and §w 2 ) and equations (A 27 ) through (A 32 ) are 
integrated (by the Adams -Moulton numerical integration method) to sufficiently 
high values of t] to insure f"(r)) -> 0 and g * ( rj ) -*• 0 for the three combina- 
tions of initial guesses (f^, g^) , (fw 2 > gwj, and (f^, g^ 2 ) • These three 
solutions give sufficient information for determining the partial derivatives 
in equations (A 3 T) and (A38) by a finite difference. With the initial guesses 
improved by Af” and Ag r the process may he repeated until the solution is 
converged upon to the desired accuracy. The requirements for convergence for 
the data in this paper were: 

|f*H - i| <0.0005 |f"H| <0.0005 

|g(co) - l| <0.0005 |g*Hi <0.0005 

This required on the average two to three cycles of the iterative process. 
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APPENDIX B 


REAL-GAS PROPERTIES 


The real -gas properties used for this report were taken from the 
literature (refs. 6 through 15 ) and were curve fitted with polynomials for use 
in the houndary -layer equations. The references were chosen primarily because 
they presented the desired properties or sufficient information to compute the 
properties for increments of temperature which were always equal to or less 
than 100 0° K for pressures ranging from 10~ 4 to 10 2 atm. This provided suf- 
ficient input data for a least-squares polynomial curve-fit method in terms 
of the desired function, enthalpy. The references are felt to indicate rea- 
sonable property variations with enthalpy and pressure for the enthalpy range 
considered here (20 < h < 18,000 Btu/lb) . 

The tabulated data were fitted in sections, generally by 7th degree poly- 
nomials, and a deviation from the fitted points within ±2 percent was main- 
tained. The equations are written in the form 

/ \ S N. V 

p = ao + ai + a 2 + ' • • • ' + a ^ (jiy 

where 

Pr/P 

Jl == { Pm*/ 

[pM'/PrPwt a w 

as the case may he. The coefficients ao through a 7 and the h/h r range 
for which they are valid are given in tables II to V. The reference values 
(p r , h r , and n r ) are given in table I. In general, the property curves are 
continuous over the complete enthalpy range but the derivatives at the points 
of connection of two sections are not. This is one reason the property deriv- 
atives were eliminated from the boundary -layer equations. The property curve 
fits for air were taken directly from reference 23 * 
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TABLE III.- PROPERTY COEFFICIENTS FOR NITROGEN 






TABLE IV.- PROPERTY COEFFICIENTS FOR HYDROGEN 



100. .01 2.21 . 2760835E 01 -.1818611E 01 -.1082485E-00 ,3076645 
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